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The reaction A^(7r, 2ti)N has been the subject of a number of experimental and 
theoretical investigations. The reaction is of interest in connection with various as- 
pects of chiral symmetry and its spontaneous breaking, including a determination of tt-tt 
scattering amplitude and the study of non-linear realization of chiral symmetry. The the- 
oretical work on the reaction A^(7r, 2'k)N may be broadly divided into two categories. In 
one category, calculations are based on a tree-level chiral Lagrangian and include meson 
and baryon resonances as explicit degrees of freedom [^]-|^. While these models capture 
the physical essence of the reaction under consideration, they are not consistent with the 
power-counting scheme of Chiral Perturbation Theory (ChPT) [^|TU|. The second category 
includes calculations which are fully within the framework of the baryon ChPT. In P the 
lowest order (0(g)) relativistic Lagrangian was employed, while in the calculations were 
performed in the next-to- leading order (0(g^)) relativistic framework. The power counting 
ambiguities of the relativistic formalism are discussed in 0. Heavy-Baryon Chiral Pertur- 
bation Theory (HBChPT) circumvents these difficulties |]rT],|T^. The most comprehensive 
study of the reaction N{tt,2tt)N at the threshold was carried out in [0, where a complete 
0{q^) calculation was performed in HBChPT. The purpose of this paper is to investigate 
the reaction N{'ir,27r)N beyond the threshold in HBChPT up to O(g^). Technically, the 
threshold and beyond the threshold calculations differ in two respects. First, as one moves 
away from threshold the number of Feynman diagrams increases significantly. Second, the 
evaluation of loop integrals becomes more involved. 

Our starting point is the low energy expansion of the vrA^ chiral Lagrangian in HBChPT 
as described in detail in Ref. [0. The effective chiral Lagrangian may be expanded as 



where the superscripts denote the chiral order in the standard power counting scheme pJ] , p!2 
In order to obtain amplitudes of 0{q^), the pion and the pion-nucleon Lagrangians may be 
truncated at O(g^) and 0{q^), respectively. The pion and the pion-nucleon, £^*|f parts 
of the Lagrangian are discussed in detail in P and |13[, respectively. In the absence of 
external fields and in the isospin symmetric limit, the pertinent pieces of the Lagrangian 
assume the following form: 

C^^^ = ^{u-u) + ^mlF'{U), (2) 

4^) = ^ {Ah{u ■ uf + Ah{u,u^){u^u,) + h{x+? 

+h{2{x+){u-u) + 2{xl)-{x-?)} . (3) 
4i = {ivV + QaS-u^ , (4) 
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Here, S'^ = i^'~i^a^^Vu is the spin matrix, U = = exp(ir-7r/F) is the SU(2) pion field with 
F being the pion decay constant, and N^{x) is the velocity-dependent nucleon field defined 
in terms of the Dirac spinor, ip{x), as 
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where M is the nucleon mass, and the four- velocity satisfies the constraint v"^ = 1. The 
last term in involving the Low Energy Constants (LECs) hi is usually referred to as the 
counterterm Lagrangian. Some of the coefficients U and hi are divergent in order to absorb 
the divergences of the one-loop diagrams which start contributing at order 0{q^). In the 
dimensional regularization scheme, these low energy constants U and hi can be decomposed 
as [0: 
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where /[(/i) and are the renormalized values of k and hi at the scale yU. Ecker and 

Mojzis [|13| have worked out the counterterms and their /5-functions involved in C^^]^ and 
we refer to their paper for the explicit form of the operators Oi and the coefficients jSi. The 
coefficients ji are given in P,p!5| . We have also used the standard notation for the operators 
X = {rriu + rridjBo = ml, = i{u^d^u - ud^u^), V = d + T and T^, = \{u^d^u + ud^u^). 



By expanding the field U 



in powers of the pion field and substituting the result in the 



effective Lagrangian, we obtain an explicit form of the Lagrangian suitable for deriving the 
Feynman rules. To a given chiral order, the transition amplitude Tfi receives contributions 
from all the possible connected Feynman diagrams up to that order. The (non-vanishing) 
topologically distinct diagrams of up to 0{q^) contributing to the reaction N{tt,2tt)N are 
given in Fig. 1. The self-energy diagrams are not shown. It is worth emphasizing that the 
actual number of diagrams for any particular physical channel is much greater than that of 
the ones displayed in Fig. 1. The individual Feynman graphs of 0{q^) are not necessarily 
finite, but the total transition amplitude is finite at each order. 

Before presenting our results we have to make a few remarks on the input parameters of 
our calculations. For the reaction under consideration up to O(g^) the following vertices and 
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LECs are needed: involving m.„ and F^r, with LECs (z = 1, ■ ■ ■ , 4), containing 
qa, containing M and ai through 05, and finally C)^n^ containing 24 LECs bi of which 
14 can, in principle, contribute to the reaction under consideration. Of the 14 LECs, 5 
combinations containing 9 LECs contribute to the vrA^ elastic scattering; these combinations 
have been recently determined |]T5[. The remaining LECs needed in our calculation are 
^4, ^11, ^12, ^13, ^14, ^17 • In this work, we treat the finite part of the unknown LECs as free 
parameters and denote them by bi. The numerical values of and bi are listed in Tables | 
and 0, respectively. Throughout the calculations the physical values of the quantities 171^, 
Fjr, gA and M are used. As discussed in [|l5l this amounts to corrections of O(g^) and higher. 
The corrections of O(g^) are explicitly included in the amplitudes of the same order. 
In HBChPT the transition amplitude for the reaction 



7r(g) + N{p = Mv + k) — > n{qi) + 7r(g2) + N{p' = Mv + lj 
can be written as 



(10) 
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in which m["^(/) denotes the heavy baryon spinor with residual momentum / and spin pro- 
jection a. We find that the amplitude A has the following Lorentz structure 
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where Ai (i = 0, 1, 2) and B are invariant functions of external momenta. 

In order to make contact with the experimental data, we need to relate the frame- 
dependent HBChPT amplitudes to the relativistic invariant amplitudes. The relativistic 
amplitude for the reaction N{tt,2tt)N can be parameterized as P,^: 

Mfi = u{p') 75 [/i + f2 4i + h 42 + h 4i 42] u{p). (13) 

The use of (|T^) and (pr3|), along with the matching condition introduced in |]TB[ results in 
the following relations between the relativistic and heavy-baryon amplitudes: 
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/i = M 1 



/2 

h 



+ 



v-qi q-q2 - v-q2q-qi + v-{q2 - qi) {ml + ^1-^2) 



--(Ao + A) - [2Mvq2 + q2-{q~qi)-ml 



--{AQ + A2) + [2Mvqi + qi-{q-q2) 



B, 



/4 = 2M 1 



4M2 



B, 



(14) 
(15) 
(16) 
(17) 



where t = {p — p') ■ At threshold the relativistic transition amplitude takes a simple form: 
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where 
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mAf2 + h)-fi-mlU ■ (19) 



The threshold amphtudes corresponding to a well-defined isospin channel, commonly referred 
to as a^'^ and a^'*^, can be decomposed in terms of physical amplitudes 0: 

ai.o = J=a-"-''"-^a-''-''^ (20) 
a3-2 = ^^-+-+n^ (21) 

where a'^"'^"" and a'^^'^^" are the threshold amplitudes for the physical reactions 7r~p vr^vr^n 
and TT^p Ti^Ti^n respectively. Furthermore, the amplitudes Di,D2 defined by Bernard, 
Kaiser, MeiBner (BKM) are related to a^'^ by 

Di = ^^'^'^ (22) 



2 a^'^ 1 

D2 = --^--a''^. (23) 

3 yiO 3 ^ ^ 



Our numerical results for Di and 1^2 up to chiral order 0{q^) are listed in Tables ^ 



and Clearly the contribution of non-leading terms to the threshold amplitudes are 
substantial, and the results seem to be sensitive to the values of the unknown LECs. Our 
calculated values of Di and D2 seem to be in fair agreement with the data. These statements 
are somewhat tentative in that the use of the actual values of the unknown LECs may lead 
to a different set of conclusions. The HBChPT calculations of BKM, on the other hand, are 
in impressive agreement with the data. In principle, the difference between our predictions 
and those of BKM is a consequence of using different values of LECs bi. In their work, BKM 
estimated the contribution of the terms containing the unknown LECs using the resonance 
saturation method As such, there are no free parameters in their calculations. It is 
conceivable that with a systematic variation of the unknown LECs we may also obtain a 
good agreement with the data. Finally, we note that the calculations of Ref . , while not 
fully consistent with the principles of chiral perturbation theory, yield results which are in 
close agreement with our calculations. 

We now turn our attention to the total cross sections. The spin-averaged cross-section 
is given by 



a 



As [{p ■ qf - M2m2]i/2 J '^/^l' (27r)3 E^, (27r)32a;,, {2'Kf2uj,, (2vr)'5'(P/ P^), (24) 



where s is a statistical factor accounting for the existence of identical particles in the final 
state. Using (0), we obtain 

[T^^ = C{p)C[p') [^^{vAvA* - \A\') + [e^p.^v-qf^qUif \B\'] , (25) 
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where C{p) = — — . To compute the cross-section, one has to integrate (p5D over an 

appropriate phase-space. The phase-space integration was carried out numerically using a 
Monte-Carlo routine GENBOD from the CERN program library |17]. It is worth noting 
that while the transition amplitude is evaluated in the isospin limit, physical masses are 
used for the reaction kinematics. 

The calculations have been performed for the Lab incident pion kinetic energy, T^, from 
threshold up to 400 MeV. From a theoretical point of view, chiral perturbation theory is 
expected to be more reliable at low energies. Hence, for each reaction, we shall present 
our results for two kinematic domains: from threshold to 200 MeV, and from threshold 
to 400 MeV. Furthermore, to facilitate comparison among different chiral orders, results 
will be displayed order by order in the chiral expansion. Finally, in order to study the 
sensitivity of the results to the values of the unknown LECs, 6j, we have performed three 
sets of calculation for each reaction: all unknown bi = 0, unknown bi = 10, and unknown 
bi = —10. The choice of the range of variation of the unknown 6j is based on the reported 
values of the known LECs bi as listed in Table ^ Some of the known LECs in Table |T| have 
rather large uncertainties. It was found that the cross- sect ions are not noticeably affected 
by these uncertainties. Our results are illustrated in Figs. 2 and 3. 

Our main results can be summarized as follows. The leading order amplitudes consis- 
tently underestimate the experimental results by a large amount over most of the kinematic 
range considered in our calculations. The inclusion of terms of 0{q^) produces only a small 
change in the cross-section. The contribution of terms of 0{q^), on the other hand, is quite 
large and in general results in an improved agreement with the data. From a phenomeno- 
logical point of view these observations can be regarded as a success for HBChPT. However, 
from a theoretical point of view, the fact that the contribution of terms of order 0{q^) is 
much larger than that of the lower order terms, is indicative of the slow rate of convergence 
of HBChPT. The slow convergence of HBChPT has also been observed 0,|l^ for processes 
other than A^(7r, 2'k)N . The lack of knowledge about the values of a number of LECs, bi, 
is the main source of uncertainty in our calculations. Our results suggest that variations in 
bi can produce quite noticeable effects in observables considered here. Clearly a systematic 
investigation is necessary in order to determine or constrain the unknown LECs. So far as 
the reaction A^(7r, 2tt)N is concerned, one could try to fit the total cross-section in conjunc- 
tion with a number of more sensitive differential cross-section and angular correlation data. 
These investigations are beyond the scope of the present work and will be the subject of a 
future publication. 

To conclude, we note that our quantitative results for the total cross-section are in fair 
agreement with previous calculations for the reaction A^(7r, 2tt)N) The fundamental 

difference between the present work and the previous ones is the calculational framework 
employed. Our calculations are performed in HBChPT and hence obey a consistent chiral 
power counting scheme. 

This work has been supported by the Natural Sciences and Engineering Research Council 
of Canada. 
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Fig. 1: Topologically distinct graphs contributing to the vrA^ — > vrvrA^ reaction up to 
chiral order 0{q^). Self-energy diagrams and diagrams which give zero contribution are not 
shown. 

Fig. 2: The total cross section in the threshold region for Tr+p n^Ti^n, n^p Tr+vr^n 
and TT^p 7r°7r°n reactions. The long dashed curve is our prediction to chiral order 0(g), the 
dashed-dotted curve corresponds to 0{q) + 0(g^), the solid curve is the complete calculation 
0{q) + 0{q^) + 0{q^) with all unknown LECs set to zero. The dotted and short dashed 
curves show the sensitivities to upward and downward shifts of all the unknown LECs by 10 
units, respectively. The experimental data are represented by the symbols: Filled diamonds: 
||T9| , empty triangles: [^, Filled squares: Filled circles: Empty circles: [^, Stars: 



H, Filled Triangles: |2| . 

Fig. 3: Same as Fig. 2 but for a wider kinematic range. 
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TABLE I. The LECs contributing to amplitudes of 0{q'^) or higher |15| 





0.2 


as 


as 


-2.60 ±0.03 


L40±0.05 


-1.00 ±0.06 


3.30 ±0.05 


TABLE IL The LECs contributing to amphtudes of 0{q^) ||,|T5| 


h 


h 


h 


u 


-2.3 ±3.7 


6.0 ± L3 


2.9 ±2.4 


4.3 ±0.9 













^1 ±62 


h 




^15 - ^16 


^19 


2.4 ±0.3 


-2.8 ±0.6 


1.4 ±0.3 


-6.1 ±0.4 


-2.4 ±0.2 



TABLE III. Calculated values of Di and D2 in this work 





0{q) 


0(g2) 


0(5'), {h = 0) 


0(g3), (6, = 10) 


0(g3), (6, = -10) 




2.50 


2.14 


2.04 


0.50 


4.46 


I?2(fm3) 


-7.70 


-6.96 


-10.88 


-8.92 


-13.95 



TABLE IV. Values of Di and D2 





Experiment yiS]] 


Present Work 
0(9'), (fti = 0) 


Theory 


Theory |Q 


Di(fm^) 


2.26 ±0.12 


2.04 


2.65 ±0.24 


1.87 


1)2 (fm^) 


-9.05 ±0.36 


-10.88 


-9.06 ± 1.05 


-10.58 
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